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ABSTRACT 

An interesting probe of the nature of dark energy is the measure of its sound speed, 
Cs- We review the significance for constraining sound speed models of dark energy 
using large neutral hydrogen (HI) surveys with the Square Kilometre Array (SKA). 
Our analysis considers the effect on the sound speed measurement that arises from 
the covariance of Cs with the dark energy density, Ode, and a time- varying equation of 
state, u;(a) — wq + {1 — a)wa- We find that the approximate degeneracy between dark 
energy parameters that arises in power spectrum observations is lifted through red- 
shift tomography of the Hl-galaxy angular power spectrum, resulting in sound speed 
constraints that are not severely degraded. The cross-correlation of the galaxy and 
the integrated Sachs- Wolfe (ISW) effect spectra contributes approximately 10 percent 
of the information that is needed to distinguish variations in the dark energy param- 
eters, and most of the discriminating signal comes from the galaxy auto-correlation 
spectrum. We also find that the sound speed constraints are weakly sensitive to the HI 
bias model. These constraints do not improve substantially for a significantly deeper 
HI survey since most of the clustering sensitivity to sound speed variations arises from 
z<1.5. A detection of models with sound speeds close to zero, Cs< 0.01, is possible 
for dark energy models with ■w> — 0.9. 

Key words: cosmological parameters - large-scale structure of the universe - cosmic 
microwave background - radio lines: galaxies 



1 INTRODUCTION 

The observed acceleration of the cosmic expansion has 
challenged our understanding of the composition and evo- 
lution of the universe. Evidence supporting the idea that 
about two-thirds of the energy in the universe is in the 
form of dark energy driving this accele ration has arisen 
from observations of ty pe la supernovae (|Riess et al.lll998l : 



tional evidence as that for dark energy, the simplest 
theoretical explanation is that of the cosmological con- 
stant (A), whose main difficulty is the dramatic incon- 
sistency between its measured v alue and the predicted 
value from quantum field theory l| Weinberg 1 1983 : ICarrollI 
[200l|). The third category includes dynamical dark energy 
theories such as quint e ssence (e.g . [Peebles fc Ratral Il988l: 
Ratra fc PeeblesI Il988l: IWetterichI Il988l: iFerreira fc Joyce 



iPerlmutter et a .1 Il999l') and from a combination of large- Il997l: ICaldwell. Dave fc Steinhardtl 11998 . also see 



scale structure ( Tegmark et al.ll2006f) and cosmic micr owave 



background (CMB) observations (|Spergel et al.|[2007l ') 



Lindei 



Theoretical explanations of the observed acceler- 
ation can broadly be classifi ed into three categories 
(|Bean. Carrol fc TroddenI |2005| ). Firstly, there are mod- 
els that remove the need for a new exotic compo- 
nent by seeking alternatives to dark energy, for ex- 
ample , by m odifying gravity o n co s mological scales 



2007^ and references therein) and k-es s ence models 
jArmendariz-Picon. Mukhanov fc Steinhardtl l200d . 1200 ll : 
IChiba. Okabe fc Yamaguchill2000l : IChiball2002l ). 



(e.g. lDvali~Gabadadze fc Po r ratil l2000l: [Pefi^ave 



Deffavet. DvaU fc Gabadadzd 1200 



I2OOOI: iDefi^aved 
^: ICarroll et al.l 



2001 



2005 



Capozziello et al.l |2006| . also see iBean et al.l (2007|) and 



references therein). Alternatively, if one assumes the va- 
lidity of general relativity and interprets the observa- 



In dynamical dark energy theories, dark energy is 
typically modelled as a scalar field, with k-essence dif- 
fering from quintessence in that it has a non-canonical 
kinetic term in the Lagrangian. The mechanism by which 
quintessence comes to dominate at later times, generating 
the accelerated expansion, is not clearly identified, leading 
to the so-called coincidence problem. The k-essence models 
were proposed to explain the late-time domination in 
a more natural way. Recently it has been claimed that 
k-essence models are not physical because in order to 
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solve the coincidence problem in these models the speed 
of sound of dark energy, c^ = 5pl5p, has to be greater 
than unity in some epoch (IBonvin. Caprini fc Durreil 
120061 . 120071 ) which violat es caus ality. Other authors, 



however, (e.g. Kang. Vanchuri n fc Winitzkil 120071 : 



iBabichev. Mukhanov fc Vikmaui i2007l ) have argued that 
superluminal sound speed propagation does not lead to 
causality violation. There also exist k-essence models in 
which the sound speed is always less than unity but these 
do not app ear to solve the coincidence problem (IScherreiJ 
l2004l20od ). 

While the viability of k-essence models is debated, it 
is nevertheless useful to consider experiments which could 
distinguish between this type of dark energy model and 
quintessence. The sound speed for quintessence is always 
unity so that, as in the case of the cosmological constant, 
the quintessence field is expected to have no significant den- 
sity fiuctuations within the causal horizon, consequently it 
should contrib ute little to the clusteri ng of matter in large- 
scale structure IjFerreira fc Jovcdl 19981 ). In k-essence models, 
however, the sound speed is not unity and dark energy can 
cluster, thereby affecting the growth of large-scale structure. 
The detection of a signature of sound speed in the integrated 
Sachs Wolfe (ISW) effect as well as in the clustering of mat- 
ter, can therefore provide valuable insight into the nature of 
dark energy. 

Many experiments have been proposed to probe the 
properties of dark energy. While many have pursued 
constraining its equation of state , w — p/ p, for example, 
using cluster counts surveys (IHaiman. Mohr fc Holderl 
I2OOII : IWeller Battve fc Knelssj* '2002"), baryon acoustic 
oscill at ions (jBlake fc Glazebrook 2003 : Hu fc Haimani 
20031: ICoora^ |2004| ) . weak lens i ng jHutereiJ I2OO2I : 
Takada fc Jaiil 12004 ISong fc K nm^ '200 4|) and t y pe la 

supernova experiments (e .g. WcUcr fc Albrechd |2002| : 

iNesseris fc Perivolaropouloj |2005| ). only a few have 
considered the prospects for sound speed detection 
llHu fc Scrantonll2004l: JPedeo. Caldwell fc Steinhardtll2004 



Weller fc LewisI 12 003*: 'Hannestadf l2005l: iBean fc Dori 
2004 ICorasaniti. Giannantonio fc Melchiorril |2005| ). In 



In this paper we extend our analysis to study the covari- 
ance of the sound speed with other dark energy parameters, 
particularly wo^Wa and fide. Since there are degeneracies in 
the parameter space when using the ISW effect and the clus- 
tering of HI galaxies as observables, it is important to con- 
sider the effect that uncertainties in these parameters have 
on the sound speed forecasts. For example, increasing the 
density of dark energy produces a larger suppression of the 
gravitational potential but this can be compensated by mak- 
ing the nature of dark energy more like matter by reducing 
its sound speed. Similarly, a lower value of w boosts the ISW 
effect at recent times but this can be compensated by reduc- 
ing the sound speed of dark energy. We implement a Fisher 
matrix analysis to study the covariance between the dark 
energy parameters and calculate the significance for ruling 
out dark energy models with constant sound speed much less 
than that of canonical models. In particular, many k-essence 
models have low sound speeds for most of the period betwee n 
last scattering and the present epoch (|Erickson et al.ll2002l ). 

We also study the sensitivity of the sound speed con- 
straints to different models for the HI bias, which has a 
direct impact on the galaxy power spectrum. Recently there 
has been discussion of radio telescopes other than the SKA 
which could carry out deep Hl-galaxy redshift surveys over 
the same area of sky, and which could potentially be built be- 
fore the multi-purpose SKA. We investigate how constraints 
on the sound speed improve for a much deeper Hl-galaxy 
redshift survey. 

The outline of the paper is as follows. In 32] we present 
the ISW and galax;y observables. In 33] we revise the proper- 
ties of our model HI redshift distribution and the HI bias. In 
33] we describe the statistical method used to derive param- 
eter constraints and its numerical implementation. Finally, 
the results are presented and discussed in [JS] with conclu- 
sions drawn in 30] 



2 THEORY 



Torres-Rodriguez fc CressI (|2007l . hereafter TRC07), we 2.1 Observables 



studied the potential for combining data from the Square 
Kilometer Array (SKAq) with forthcoming Planclo data to 
measure the sound speed. Assuming the SKA will provide a 
redshift survey of 21-cm emitting galaxies over most of the 
sky out to z ~ 2, and a specific model for the evolution of 
HI in galaxies, we considered a combination of the galaxy 
auto-correlation and the galajcy-CMB cross-correlation in 
several cosmic epochs to trace the evolution of clustering. 
This evolution is somewhat sensitive to the sound speed 
of dark energy in the low-sound-speed regime. The cross- 
correlation signals are dominated by contributions from the 
ISW effect and thus provide a probe of the evolution of 
potentials on large scales where clustering of dark energy 
could be most noticeable. We found that for models with 
constant 10, given high precision measurements of other 
cosmological parameters, it was possible to measure the 
effects of a sound speed that is much smaller than unity. 



^ http://www.skatelescope.org/ 

^ http : //www . r ssd . esa . int /index . php?pro j ect=PLANCK 



We first review the ISW and Hl-galaxy clustering observ- 
ables which are used to discriminate dark energy models. 
More details can be found in TRC07. 

The fiuctuations in the matter distribution of the large- 
scale structure can be expressed in terms of the projected 
fractional source count of the mass tracer 

-(n) = / bHi[z) -7—dm[z,n)dz 



No 



dz 



dN . 



5m{n) I bHi{z)-—D{z)dz, 
Jo dz 



(1) 



where bni is the linear bias parameter of the Hl-galaxy pop- 
ulation, 5m is the matter overdensity {5^ = 5m{z = 0)), 
dN /dz is the normalised redshift distribution of HI galaxies 
and D{z) is the linear growth of matter fiuctuations given 
hyD{z)=5m{z)/5'L. 

The Fourier modes of the temperature fiuctuations orig- 
inating from the ISW effect are expressed as the change in 
the gravitational potential over conformal time (or comov- 
ing distance, r) integrated from today to the epoch of de- 
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coupling 



ST 



(fc) 



dr <l>'(r-, fc) 



e2fc2 ^"^(fe) 



rfff(g) 

dz 



dz . 



(2) 



where _ffo and Sim are, respectively, the value of the Hub- 
ble constant and the matter density parameter today, $ is 
the Newtonian gravitational potential, 5%^{k) is the Fourier 
transform of the matter distribution and the prime denotes 
derivatives with respect to comoving distance. The domi- 
nant contribution to the ISW effect comes from the CDM 
perturbations. This allows us to express the evolution of the 
gravitational potential as the change of the linear growth 
suppression factor, g{z) = (H- z)D[z), via the Poisson equa- 
tion. 

We can express the galaxy auto-correlation, CMB auto- 
correlation and galaxy-CMB cross-correlation in harmonic 
space via their respective angular power spectra as 



Cf 



c 



c. 



gT 



^ , dk /i5r,,,|2\ ,. 

= 4^ / ^/f(fc)/7(fc)A^(fc), 



(3) 



(4) 



(5) 



where ji{kr) is the spherical Bessel function and Am(fc) = 
k^Ps{k) /2tv^ is the logarithmic matter power spectrum to- 
daywithPa(fc) = (|5i(fc)p). The functions /i^(fc) and/f (fc) 
correspond to the weight functions for the HI survey and the 
ISW spectra respectively. They depend on both the redshift 
distribution of the galaxy selection function and the rate of 
change of the gravitational potential according to the cos- 
mological model of dark energy. From Eqs. ([l][2} the weight 
functions are defined as 



/f(fc) = 



bHi{z)-j-D{z) jt{kr{z)) dz , 



dz 






dz 



ji{kr{z))dz. 



(6) 
(7) 



2.2 The effect of dark energy on the observables 

The background expansion of the Universe is affected by 
the density of the matter component and by the equation of 
state parameter, w, through the Hubble parameter, H{z). 
This modifies the comoving distance, r, and has a direct 
effect on Eqs. dMSll. 

In a general fluid description for the density perturba- 
tions of dark energy (e.g. Bean & Dore 2 0041 ') , the evolution 
of the fluctuations is characterised by both the equation of 
state and the speed of sound. In the frame where cold dark 
matter (CDM) is at rest, the evolution of the density and ve- 
locity perturbation of a general matter component (denoted 



by subscript i) is given by (|Weller fc Lewisll2003h 

5- -f- 3W(c^,, - w,){5, + 3?i(l + w,)vi/k) + 

(1 + Wi)kvi + SHw'v/k — —3(1 + Wi)h' 

v[ + H{1 - 3cs,i)wi = fcc^,i5i/(l + Wi) , 



(8) 
(9) 



where Ti is the conformal Hubble parameter, Vi is 
the velocity, h' the sy nchronous metric perturbation 
(IMa fc Bertschingerlll995l ) and the circumflex (") means the 
sound speed is deflned in the rest frame of the dark energy 
component. Eq. ^ includes the variations of w with respect 
to conformal time. We consider models of dark ener gy with 
a slo wly varying equation of state parameterized IjLindeii 
[2003) as a function of the scale factor, a: 



w{a) = Wo + [1 — a)wa. 



(10) 



We note that ildc only appears directly in the back- 
ground expansion, unlike w(a) which, together with Cs, has 
a direct contribution to the density evolution. As we will see 
later, this has an effect on the covariance of these parame- 
ters. 



3 THE HI-GALAXY SURVEY 

In TRCQ7 we presented the motivation for using an exper- 
iment like the SKA to complete a large HI survey. There 
are some constraints on the distributio n of HI gas at high 
redshift from Ly-a absorption studies (jPeroux et al.l 12003! ) 
but there remains a great deal of uncertainty in the HI se- 
lection function at high redshift. As our refer e nce m odel, we 
adopt the model 'C in lAbdalla fc Rawlingsl (|2005l ) and in- 
vestigate the impact of changing the selection function by, 
firstly, changing the integration time and, secondly, changing 
our assumption about the evolution of the HI bias. 

For a single-pointing survey, a field of view (FOV) fre- 
quency dependence of i/~^ and a detection threshold of 
S/N — 10, the number of HI galaxies per square degree, 
per redshift interval and at redshift z is expected to be rea- 
sonably approximated by 



dN 

— ;— — Az exp 

dz 



2 a? 



(11) 



where the fitting parameters A, Zc and a^ depend on the in- 
tegration time. We would like to investigate the performance 
of two different surveys: one of 4 hours of integration time 
per pointing (as in TRC07) and an ultra deep survey, using 
36 hours of integration. We assume that both surveys are 
completed by an SKA-like experiment and cover the same 
area of sky but that the latter survey takes nine times longer 
to complete. Both survey selection functions are depicted in 
Fig. [11 

The selection function of our mass tracer is then 
0(2;) = bni{z) X dN /dz, where the tilde represents the nor- 
malised distribution over the redshift of interest. In order 
to increase the significance of distinguishing models, we di- 
vide the galaxy selection function into several redshift bins 
of width Az = 0.2 up to Zmax = 2 (or Zmax = 3 for the 
deeper survey), the width of redshift bin being chosen to 
minimise shot noise from a small number of galaxies in nar- 
row bins and smearing of the gravitational potential along 
the line of sight in wide bins. We consider the galaxy and 
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SKA 4hrs 
SKA-like 36hrs 




Figure 1. SKA redshift distributions of HI galaxies and models 
for the HI bias parameter used. An integration time for a single 
pointing of 36 ho urs (black dashed line) would reach a depth of 
2 = 3. Taken from lBasilakos et al.l 1120071) . the two models for the 
HI bias (in grey) only differ greatly at high redshift. A constant 
bias model (grey dashed line) is unlikely to be realistic. 



the ISW-galaxy power spectra in each bin as independent 
measurements. 

From Eq. ((6)l, we see that the bias factor has a direct 
impact on the galaxy-galaxy spectra. The effect of the bias 
model on the sound speed constraints is not straightforward 
to see, so we investigated this further by using different bias 
models. We have used two mod els for bni based on the 
studies of iBasilakos et al.l ((20071) (and references therein). 
In the first model we make the simple but unrealistic as- 
sumption that the HI bias does not evolve with redshift and 
that it r emains constant at the value measured locally by 
HIPASS (IBarnes et al. 2001 ) in th e concordance model, i.e., 
hni = 0.68 ( Basilakos et al.ll2007l l. More realistically, a sec- 
ond HI bias model considers the transformation of gas into 
stars and the gas consumption as a function of redshift. In 
this model HI galaxies are more biased at high redshift rel- 
ativ e to the local va l ue. U nder the linear perturbation the- 
ory, IBasilakos et al.l (|2007l ) derive the HI bias evolution to 
high redshift for a ACDM model, which is shown in Fig. [T] 
together with the constant bias model. We note that an in- 
dependent measure of the bias of Hl-selected galaxies as a 
function of redshift will be possible using redshift-space dis- 
tortions in the survey. 



4 STATISTICAL METHOD 

4.1 Numerical implementation 

We compute the relevant power spectra in Eqs. (|3H5J us- 
ing a modified version of the CROSS_CMBFASt|j code. These 
modifications include the addition of our Hl-galaxy redshift 
distribution (see ^ as well as the bias factor. The result- 
ing selection function is split into different redshift bins for 
which the relevant power spectra are calculated. 

We choose fiducial models that give an angular di- 
ameter distance to recombination fixed to the value from 



|http : //www . astro . Columbia . edu/~plerste/ISWcode . html | 



WMAP observations l|Spergel et al.ll2007h and let the dark 
energy parameters f2de, wo, Wa and c^ vary around the fidu- 
cial model. Unless otherwise indicated, the dark energy fidu- 
cial model is: fido = 0.705, wo = —0.8, uia = and Ca = 1. In 
addition to the dark energy parameters, we assume a phys- 
ical matter density of Q.mh^ = 0.126, a physical baryonic 
density of Q.bh^ — 0.0223, an optical depth to reioniza- 
tion of r = 0.09, a primordial power spectrum amplitude of 
A|> = 2.02 X 10~^ a nd a scalar spectral index of n^ = 0.951 
(|SpergeI et al.ll2007l ). 



4.2 Forecasts for sound speed detection 

In order to measure the significance of a detection of a sound 
speed model we utilise the Fisher information matrix. The 
Fisher matrix is defined as 

where C, the likelihood, is the probability of observing the 
data set {xi,X2, ■■■} for a given cosmological parameter set 

{^1,^2,...}. 

The Fisher matrix method allows us to forecast how 
well a survey will perform in constraining a set of cosmo- 
logical parameters, by providing the minimum systematic 
uncertainty in each model parameter that is to be fit by the 
future survey data. Under the assumption that the individ- 
ual parameter likelihoods approximate a Gaussian distribu- 
tion, the informati on contained in the angular power s pectra 
can be written as IjTegmark. Tavlor fc Heavenslll997l ) 

(21 + 1) 



-Fq/3 — /sky 2_^ ■ 



■ Tr[DtoC7"D£,3C7 



(13) 



where the sum extends over multipoles £, fsky is the amount 
of sky covered by the survey, Ce is the data covariance ma- 
trix, and T)ea is the matrix of derivatives of the angular 
power spectra 

D--??| (14) 

with respect to the parameters, 6a , evaluated at the fiducial 
model. The data covariance matrix elements include the an- 
gular power spectra plus the noise terms. For the CMB, the 
noise contribution to the temperature measurement depends 
on the beam window function and the pixel noise of the ex- 
periment 






lp-2 
T -"« 1 



(15) 

where Bp is the window function of the Gaussian beam 
and Wy^ = (Tp ^boam is the inverse noise weight with 6'bcam 
and cTp, respectively, the beamwidth and noise per pixel in 
a given frequency band. We consider a Planck- like CMB 
experiment that measures temperature anisotropies in two 
high frequency bands, 143 and 217 GHz. The d etails of these 
parameters for Planck can be found in e.g. iRocha et al.l 
(2004) and are not presented here. As shown in TRC07, 
the contribution from the CMB spectrum to the detec- 
tion significance is much less important compared to the 
cross-correlation and galaxy auto-correlation spectra. For 
the galaxy field, the source of noise comes from Poisson fluc- 
tuations in the number density 



Cf 



Cr + l/ni , 



(16) 
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where n^ is the galaxy number per steradian in the redshift 
bin of interest. 

We construct our data covariance matrix by combining 
the observables as follows 









c. 



gT,2 



\ Ci 



gT.lO 
















... cf^- \ 







cr''" 



(17) 



where the superscript number refers to the redshift bin and 
the zeros indicate no cross-correlation between the galaxy 
power spectra from different redshift bins. We have also as- 
sumed that there is no correlated noise between the galaxy 
and CMB power spectra. 

Information about constraints on the cosmological pa- 
rameters is contained in the derivatives in Eq. (I13|l . We cal- 
culate these derivatives for each multipole by fitting a curve 
through the power spectra values as they vary as a func- 
tion of the dark energy parameters wq, Wa and Qdc- We 
then numerically calculate the slope at the fiducial point. 
We require a precision that is less than one percent in order 
to measure the small variations in the cross-correlation and 
auto-correlation spectra. For this reason, we look at changes 
in the parameters that are large enough to allow a polyno- 
mial curve fitting. 

We present the derivatives of the auto-correlation and 
cross-correlation power spectra with respect to the dark en- 
ergy parameters in Fig.[2]and Fig.|3] respectively, for a series 
of redshift bins. We have multiplied the derivatives by the 
factor {2£ + 1)^ " to make the contribution of this product 
to the Fisher information matrix as a function of £ more 
evident. We note that the galaxy auto-correlation spectra 
contribute most of the signal in distinguishing sound speed 
models. This is due to the larger amplitude of the galaxy 
auto-correlation derivatives with respect to the dark energy 
parameters as compared to the cross-correlation derivatives. 
As we will see in the next section the cross-correlation spec- 
tra only add about ten percent of the discriminating signal. 
The cross-correlation signal originates from the large-scale 
ISW effect, thus the information from the cross-correlation 
spectrum is mainly available at low multipoles, £ <^ 20, as 
can be seen in Fig. O 

The power spectrum derivatives also provide insight 
into the covariance between the sound speed and other dark 
energy parameters. We observe in Fig. [2] that the deriva- 
tive with respect to f2dc is two orders of magnitude larger 
than the derivative with respect to the sound speed. The 
larger amplitude and distinct shapes of the ildc derivatives 
compared to the sound speed derivatives in different red- 
shift bins breaks the degeneracy between fide and c^. The 
amplitude of the change in the galaxy auto-correlation and 
cross-correlation spectra due to a varying equation of state 
(wo and Wa) is smaller than the change due to the dark 
energy density but still an order of magnitude larger than 
that caused by variations in c^. In addition the characteristic 
shapes of the auto-correlation and cross-correlation deriva- 





100 1000 



-5x10 




10 



100 1000 



/ 



Figure 2. Derivatives of the galaxy auto-correlation power spec- 
tra with respect to the dark energy parameters Hdc, u>a, wq and 
c2 in redshift bins 0.0 < z < 0.2 (solid), OA < z < 0.6 (dashed), 
1.0 < 2 < 1.2 (dot-dashed) and 1.4 < 2 < 1.6 (dotted). The 
fiducial model used is wq = —0.8, Wa = 0, Ode = 0.705 and 




Figure 3. Derivatives of the galaxy-ISW cross-correlation power 
spectra with respect to the dark energy parameters fidci ^a, wo 
and c^ in redshift bins 0.0 < 2 < 0.2 (soUd), 0.4 < 2 < 0.6 
(dashed), 1.0 < 2 < 1.2 (dot-dashed) and 1.4 < 2 < 1.6 (dotted). 
The fiducial model is as in Fig. [2] 



fives with respect to wo and Wa in different redshift bins 
allow the degeneracy between wo,Wa and Cg to be lifted. 

We extend the sum in Eq. (|13p up to ^max = 1500 to 
include the galaxy auto-correlation signal from the highest 
redshift bins, which is projected onto small angular scales. 



4.3 A quantitative approach 

We have seen how the use of a survey split in redshift bins 
helps lift the degeneracies that exist between dark energy 
parameters. We wish to quantify the effect on the sound 
speed measurement by marginalising over Wa, wq and fide- 
The full Fisher matrix for this experiment includes matrix 
elements for all parameters 
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5 RESULTS 



100 



no marginalisation 
varying a^ 
varying w(a) 
full marginalisation 




-0.85 -O.i 

fiducial w(a) 

Figure 4. The significance of a detection of a c^ = dark energy 
model compared to a c^ = 1 quintessence model as a function 
of the equation of state, w{a). The S/N for a 4-hour integration 
per pointing is shown with all parameters fixed to fiducial values 
(black dashed), varying 0^^ only (grey solid), varying w(a) only 
(grey dashed), and full marginalisation (black solid). 



i^ cc 



B 
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(18) 



where 

A = 

and 
B = 



( F„ 



Fw„c -f 






-'"'"ondo 






\ 
) 



Fn^.c ) 



(19) 



(20) 



We obtain the marginalised Fisher matrix by taking the 
inverse of F and extracting the sub-matrix corresponding to 
the parameters we want to measure. For the sound speed 
this is 

(F-^)n = (F,,-BA-^B^)^\ (21) 

and the marginalised Fisher matrix is the inverse of this 



BA^B^, 



(22) 



which is a scalar, since we are interested in measuring the 
sound speed parameter. 

In order to compute the derivatives for the Fisher ma- 
trix elements involving the sound speed in Eq. (|13|l , we note 
that the effect on the angular power spectra due to changes 
in the sound speed is only significant for order-of -magnitude 
variations of Cs . Following Hu fc ScrantonI (|2004l ) , we there- 
fore define the derivative with respect to c^ in Eq. (|14p as a 
finite difference 



D,,2->C,(c^ = l)-C,(cJ/l), 



(23) 

noting that, with this definition, the value of -Fmarg can be 
interpreted as the significance, (S/N)^, of a detection of a 
sound speed model with cj 7^ 1 relative to a quintessence 
model with c? = 1. 



5.1 Signal to noise of dark energy sound speed 
detection 

The basic detection level of a c^ = model relative to a 
Cs = 1 model is shown in Fig. |4] as a function of the equa- 
tion of state, Wa- It is clear that the significance of a sound 
speed detection is larger for models with a larger constant 
equation of state. This results from the fact that a compo- 
nent with a larger constant equation of state is more like 
matter, hence variations in its sound speed produce larger 
observable effects. For models with mq > — 0.9 the squared 
signal to noise exceeds 10 suggesting that these models can 
produce a detectable effect. 

Fig. |4] also illustrates the effect of the covariance be- 
tween the sound speed and other dark energy parameters. 
We note that the full covariance degrades the sensitivity to 
the sound speed but not significantly; in the absence of co- 
variance with all other dark energy parameters the sound 
speed constraints improve by at most twenty percent. The 
lack of severe degradation in the sensitivity, that one would 
expect to arise from degeneracies between the dark energy 
parameters, is the result of information gained from measur- 
ing the auto-correlation and cross-correlation spectra in sev- 
eral redshift bins. We note that the sound speed covariance 
with the equation of state parameters, wo and Wa, degrades 
the sensitivity more than the covariance with the density 
Side, a result that was anticipated from studying the power 
spectrum derivatives in the previous section. 

The significance of detection of dark energy models with 
Cj < 1 relative to a quintessence model with c^ = 1 is shown 
in Fig. [S] for a fiducial equation of state, w{a) = —0.8. It 
is clear that the sound speed detection is most significant 
for dark energy models with Cs ^ as the dark energy 
clustering is most pronounced in these models, and becomes 
marginal for models with Cs>0.1. We also note that the 
effect of the covariance with other dark energy parameters 
on the sound speed measurement is most evident for models 
that have a high significance of sound speed detection, and 
hardly noticeable for models in which the sound speed is 
poorly measured. 

We next consider the impact of the uncertainty in our 
HI bias model on the sound speed constraints. In Fig. [6] we 
compare the detection significance as a function of the sound 
speed for our two HI bias models. The constant bias model, 
which is unlikely to be realistic, predicts lower source counts 
at high redshift. Nevertheless this model has a sensitivity to 
the sound speed that is only about ten percent worse than 
the sensitivity of the evolving bias model. This indicates that 
the uncertainty in the HI bias does not change our forecasts 
significantly. 

Finally, we have explored how the survey depth affects 
the significance of sound speed detection. In Fig. |6] we con- 
sider the significance of detection for a 4-hour-per-pointing 
survey compared to a 36-hour-per-pointing survey. We have 
assumed that both surveys cover the same fraction of the 
sky so that the 36-hour-per-pointing survey takes nine times 
longer to complete. The longer survey is unlikely to be prac- 
tical in terms of total integration time but provides a useful 
guide as to how significant a much deeper survey will be for 
constraining the sound speed. From Fig. |6] it is clear that 
the deeper survey is able to discriminate more easily be- 
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Figure 5. The significance of separation between dark energy 
models with c^ ^ 1 and quintessence (c^ = 1) for a fiducial 
w = —0.8 model of dark energy. The S/N for a 4-hour integration 
per pointing is shown with all parameters fixed to fiducial values 
(black dashed), varying 0^^ only (grey solid), varying w(a) only 
(grey dashed), and full marginalisation (black solid). 



tween sound speed models. The signal to noise increases by 
a factor of two to three for low sound speed models but is 
approximately the same for models with Cs — > 1. 

It is interesting to ask whether this improvement arises 
from the measurement of clustering at higher redshifts, 
z>2, or from the increased number counts at intermedi- 
ate redshift. In Fig. [7] we plot the cumulative contribution 
of different redshift bins to the discriminating signal for a 
Cs = and w = —0.8 model. We note that for the fiducial 
4-hour-per-pointing survey nearly all the signal accumulates 
hj z = 1.5, and only ten percent of the signal comes from 
the cross-correlation spectra. For the deeper survey, approxi- 
mately eighty percent of the information comes from z <J 1.5. 
This suggests that the improved signal in Fig.[6]arises mostly 
from the increase in the number counts at intermediate red- 
shifts (from z — 0.5 to z = 1.5), which results in a lower vari- 
ance in Eq. 1131 rather than the clustering signal at higher 
redshift. 



6 CONCLUSIONS 

We have studied the potential of large HI surveys to con- 
strain constant sound speed models of dark energy. We in- 
vestigated the covariance between the dark energy cosmo- 
logical parameters, finding that uncertainties in the density 
of dark energy and in its equation of state will not dramat- 
ically degrade our ability to detect the sound speed. This 
arises because of the ability of these surveys to detect large 
numbers of HI galaxies in several redshift bins. The slight 
reduction in the signal to noise comes mostly from variations 
in the equation of state parameters. 

We have also investigated the impact of using an ul- 
tra deep SKA-like HI redshift survey and assessed the effect 
of changing our HI bias model. We discovered that a deep 
survey of HI galaxies, with 36 hours of integration time per 
pointing improves constraints on the sound speed as com- 
pared with a smaller 4-hour-per-pointing survey due to the 



Figure 6. The significance of separation between dark energy 
models with c^ ^ 1 and quintessence (c^ = 1) for a fiducial 
w = —0.8 model of dark energy. The marginalised S/N for an 
SKA 4-hour-per-pointing survey (solid line) and a deeper sur- 
vey of 36 hours per pointing integration time (dotted line). Both 
curves above are for an evolving HI bias model. The dashed line 
corresponds to a constant HI bias model. 



2.5 



1.5 



0.5 



total signal (4 tir) 

galaxy autocorrelation (4 hr) 

. — total signal (36 hr) 

. . . ■ galaxy autocorrelation (36 tir) 




Cs=0 

fiducial w(a) = -0.8 
2 2.5 



Figure 7. The cumulative significance S/N as a function of 
redshift contributing to the sum in Eq. II13I I for a c^ = and 
w = —0.8 dark energy model. The curves are normalized to the 
maximum value of the total signal from the 4-hour-per-pointing 
survey. The dashed line represents the signal from the galaxy 
auto-correlation alone while the solid line represents the total sig- 
nal from the auto-correlation and cross-correlation spectra. Both 
curves arc for the 4-hour-per-pointing survey. The corresponding 
curves for the deeper 36-hour-per-pointing survey are shown as 
a dotted line (auto-correlation alone) and as a dot-dashed line 
(both auto- and cross-correlation). 



increased number counts at z<,1.5. A maximum redshift 
depth of Zmax ~ 1.5 provides most of the discriminating sig- 
nal for both surveys. In addition, our results have not shown 
a strong dependence on the uncertainty in the HI bias model. 
These results could guide future planning for these types of 
survey experiments. 

Regarding the detection of the sound speed, we found 
that we can only detect models of dark energy with small 
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values of the constant sound speed, Cs <0.01. As c^ — + 0, 
a model with w — —0.9 can be detected at the 3-cr level. 
For larger values of wo sound speeds closer to zero can be 
detected with greater confidence. The study of the cluster- 
ing properties of dark energy through its sound speed thus 
promises to be an interesting approach to confront the pre- 
dictions of theoretical dark energy models and uncover the 
nature of this mysterious component. 
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